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摘 要 :研究 群 组 中 个 体 之 间 的 非 定 常 流动 机 理 ,可 以 为 仿生 飞行 器 集群 运动 提供 理论 基础 。 采 用 

de 流体 力学 方法 ,对 前 飞 状态 的 扑 翼 群 组 个 体 之 间 水 平 间 距 对 气动 性 能 的 影响 进 
行 研究 。 研 究 发 现 ,水 平 间距 对 扑 翼 气 动 性 能 具有 显著 影响 。 在 一 定 的 垂直 间距 下 ,和 群 组 中 扑 翼 可 

:之 以 在 较 小 的 水 平 间 距 下 获得 最 佳 的 推力 性 能 ,在 较 大 的 水 平 间距 下 可 以 获得 最 佳 的 升力 性 能 。 扑 

刁 翼 气动 性 能 的 变化 主要 与 群 组 中 前 翼 和 后 翼 的 脱落 涡 相 互 作用 密切 相关 。 
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Effect of individual horizontal spacing on 
aerodynamic performance of flapping wings in flocking 
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Abstract;The study of the unsteady flow mechanism between the individuals in a group can provide theo- 
retical basis for the swarm motion of bionic aircraft. The influence of horizontal spacing between flapping 
wing flockings on aerodynamic performance in forward flight is studied by using computational fluid dy- 
namics method based on finite element method. It is found that the horizontal spacing significantly affects 
the aerodynamic performance of flapping wings. Under a certain vertical spacing , flapping wings in flocking 
can obtain the best thrust performance at a small horizontal spacing, and the best lift performance can be 
obtained at a large horizontal spacing. The aerodynamic performance of flapping wings is closely related to 
the vortex shedding interaction between the front wing and the rear wing in flocking. 


Key words; aerodynamic performance ; flow mechanism; multi-flapping wing; LEV 
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在 微型 飞行 器 等 仿生 技术 的 推动 下 ,生物 飞行 
力学 研究 受到 众多 学 者 的 关注 , 人 们 已 经 对 单个 扑 
BAT TRAIL ,但 在 自然 界 中 生物 集群 现象 
是 普遍 存在 的 ”1 。 进 一 步 研究 群 组 之 间 个 体 非 定 
常 流动 机 理 , 将 有 助 于 促进 仿生 工程 的 发 展 。 

部 分 理论 和 实验 研究 "证明 : 鱼 群 和 鸟 群 通 
过 采用 不 同 的 行为 策略 ,可 以 有 效 地 利用 涡 来 提高 
其 动力 学 性 能 , 从 而 减少 能 耗 。 学 者 通过 采用 规则 
排列 的 多 个 拍 动 翼 来 评估 群 组 动力 学 性 能 。 包 含 2 
个 拍 动 翼 的 串联 和 并 联结 构 , 已 经 得 到 广泛 的 研究 。 
研究 发 现 , 拍 动 翼 间 距 、 相 位 差 和 斯 特 劳 哈 尔 数 ( Sr) 
对 群 组 中 个 体 的 推力 性 能 有 较 大 影响 5*21 。 例 如 ， 
在 串联 排 布 的 两 个 拍 动 翼 中 , Boschitsch 等 发现 
下 游 拍 动 受 由 于 受到 上 游 拍 动 慢 尾 流 的 影响 ,推力 
和 推进 效率 可 达 单 翼 的 1.5 倍 或 0.5 倍 ,其 动力 学 
性 能 与 间隔 间距 和 相位 差 相 关 。Deng 等 ”发 现下 
Die SBE LEBER Sr 下 会 使 推力 得 到 增 大 ,而 在 相对 
FES Sr 下 则 相反 。 同 时 在 两 个 拍 动 翼 并 联 时 , 同 相 
拍 劲 , 拍 动 翼 以 减少 推力 为 代价 获得 推进 效率 的 提 
项 而 反 相 拍 动 , 拍 动 辟 则 表现 为 推力 增强 ,但 推进 
效 雍 没有 明显 变化 ”” 。 
而 自然 生物 集群 中 运动 构 型 更 加 多 样 。 研 究 表 
时 不 同 运 动 构 型 ,个体 间距 和 拍 动 频率 对 个 体 和 和 群 
组 的 动力 学 性 能 都 有 显著 影响 2 。Han SE") 研 
完了 种 不 同 构 型 的 三 维 多 翼 推 力 性 能 发 现 :在 三 角 
和 过 形 结构 中 ,个 体 间距 越 近 , 群 组 整体 的 推力 表现 
越 妹 。 同 时 关于 鸟 类 的 集群 运动 已 经 有 学 者 通过 理 
论 弥 析 和 生物 实验 观测 证 明 : 对 于 以 Y 形 结构 飞行 
的 镶 类 , 鸟 群 个 体 之 间 积 极 的 相互 作用 可 以 使 个 体 
最 大 化 地 捕 提 有 益 的 上 升 气流 :2 ,进而 实现 节 
能 。 

综 上 所 述 ,目前 对 生物 集群 现象 的 流体 研究 ,大 
多 数 以 研究 鱼 群 的 推力 性 能 为 主 ,采取 的 运动 形式 
被 简化 为 翼 型 的 沉浮 运动 ;而 对 于 采用 昆虫 或 鸟 类 
运动 方式 的 集群 现象 ,同时 考虑 推力 和 升力 的 研究 
目前 还 很 少 , 且 对 于 鸟 类 集群 飞行 的 研究 多 采用 理论 
分 析 和 生物 实验 观测 ,无 法 给 出 具体 的 流 场 信息 。 本 
文 研究 个 体 水 平 间距 对 群 组 中 扑 费 气动 性 能 的 影响 ， 
并 从 扑 细 脱落 涡 的 角度 ,分 析 了 集群 运动 的 流动 机 理 。 


1 扑 曼 模型 和 运动 方式 


扑 滨 简化 为 矩形 平板 ,其 厚度 为 3% GHAR (e), 
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RIZE (R/c; R FET RRR) , 疲 根 位 于 矩形 平板 的 
一 角 ” ,如 图 1(a) 所 示 。 图 1(b) 和 图 1(c) 描 述 了 
扑 到 的 运动 方式 。 扑 既 的 运动 由 两 部 分 组 成 : 拍 动 
(HE X Ht; ULE 1 Cb) ) 和 旋转 ( 绕 扑 翼 前 缘 "” ) 。 
其 中 扑 翼 拍 动 由 下 面 的 方程 定义 , 即 

b(t) =P, +0.5@,sin(2a/t) (1) 
式 中 :g(t) 是 拍 动 角 ( 前 缘 和 Z 轴 之 间 的 角度 ,图 1 
(b) ) sf 是 扑 嗓 拍 动 频率 ;| 是 拍 动 振幅 ; ,是 平均 
拍 动 角 。 扑 驾 的 俯仰 角 运 动 方程 (a, 既 将 与 X HZ 
间 的 角度 ) 为 

a(t) =au+aisin(2T 丰 +D) (2) 
SUP: a, ETA Tie HS 5 an EF fA sp ETA AT 
(iz cya AE A HP , =90° 时 已 被 证 明 可 
以 产生 最 佳 推 力 性 能 。 


前 缘 


(b) 扑 翼 运动 示意 医 


0 0.25 0.50 0.75 1.00 
(© saz 12 cp AUF cy NP BEHET 
到 1 扑 避 模型 及 运动 方式 

Fig.1 Flapping wing model and mode of motion 

TEVA Pi HD, Dn a 和 a 分 别 为 100°、 
30° 20" 和 15°, Ð, P WRU SAE WERT ESE 
验 中 的 测量 数据 相似 。 选 择 w 和 a, 的 值 是 为 了 
保证 产生 正 升力 和 推力 ,这 与 Broering 和 Lian”! 使 


用 的 值 相似 。 来 流速 度 (V) 沿 着 X 轴 正 向 ( 见 图 
2)。 推 进 比 (= V2BR) 为 0.3, 相 当 于 真实 昆虫 
的 中 等 前 进 速度 


到 2 扑 惨 的 集群 编队 构 型 示意 图 
Fig.2 Configuration diagram of flapping wing group formation 
群 组 中 所 有 扑 慢 采 用 相同 的 运动 ,同时 定义 了 
2 个 参数 (HH 和 7y) 来 描述 群 组 中 扑 咽 的 相对 位 置 。 
如 图 2 所 示 , 水 平 间 距 及 是 相 邻 2 PS Fh Be TAR a 
的 人 间距 ,垂直 间距 y 是 来 流 方向 与 3 HS FDR ( $8 
动 赴 心 ) 连 线 之 间 的 夹 角 。 扑 翼 群 组 中 3 个 扑 翼 从 
加 依次 定义 为 Wl 、W2 和 W3。 


2 忆 计 算 方 法 及 验证 


控制 方程 和 数值 方法 


€ EIFE D EEA n] Eä Navier-Stokes 方程 
eR tA 一 化 形式 为 

V:-u=0 (3) 
Èl tu V iuz -Vp+ (Rv (4) 


R 
Nai : U 是 量 纲 为 一 的 流体 速度 ;p 是 量 纲 为 一 的 流 
Ee E 量 纲 为 一 的 时 间 ; VEEE V E 
HPI AE; Re 是 雷诺 数 。 采 用 有 限 元 法 对 式 
(3) ~ 式 (4) 进行 了 数值 求解 。Galerkin 最 小 二 乘 
(GLS) 公 式 为 所 有 变量 的 空间 离散 化 提供 二 阶 
精度 。 采 用 任意 拉 格 朗 日 - 欧 拉 (AILE) 09 方法 处 理 
RISA Nieman HK. 
KIER) rE SOON EA PR LAY 7, Fb He 

推力 (7) 定 义 为 与 来 流 平行 且 反 向 的 力 。 气 动 功率 
Pe ATI ES Fb EIS BH A REY A FD HE 
力 和 气动 功率 系数 分 别 定义 为 C, = 2L/pU’S, Cy = 
2T/pU’ S,Cp, =2P,/pUS, EP: p 是 流体 密度 ;5 
EFM Ta AR, U 是 扑 翼 回 转 半径 处 的 最 大 拍 动 速度 。 
压力 系数 定义 为 C, =2(p-p**)/pU", HP :p 是 
压力 ;p “是 无 穷 远 处 的 压力 。 升 力 和 推力 效率 分 
别 定 义 为 NL = C,/ Cp Fin = Cr/ Cpao Ci\Cr 和 Cr。 
分 别 表示 一 个 周期 内 的 平均 C, .Cy 和 Cp ao 
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2.2 求解 器 验证 和 网 格 无 关 性 验证 


为 了 验证 求解 如 的 准确 性 和 可 靠 性 ,利用 文献 
[37] 实验 测量 结果 和 文献 [38 ] 数值 计算 结果 对 本 
研究 使 用 的 求解 器 进行 验证 (注意 ,实验 测量 结果 仅 
给 出 升力 系数 ) 。 采 用 典型 的 果 蝇 翅膀 形状 ,计算 
Re =136 WIRES RS FAEERE Thaw 
采用 俯仰 角 a 提前 翻转 的 运动 形式 (具体 运动 方程 
见 文献 [38] ) 。 在 图 3 中 给 出 了 本 研究 所 用 求解 器 
数值 计算 结果 与 其 他 人 研究 结果 的 对 比 ,发现 计算 
结果 吻合 较 好 ,因此 可 以 认为 该 求解 器 可 靠 。 


当前 计算 经 


HR 
3h o SRSA SE 
o 文献 [37]Han 川 量 结 


0 0.25 0.50 0.75 1.00 
t 
(a) 


一 一 当前 计算 结果 
文献 [38] 数 值 计算 结果 


0 0.25 0.50 0.75 1.00 


图 3 不 同 数值 计算 方法 和 实验 结果 的 比较 
Fig.3 Comparison between the results calculated with 
different numerical methods and experiment 
在 计算 中 采用 非 结 构 的 四 面体 网 格 。 为 确定 合 
适 的 网 格 ,进行 了 网 格 无 关 性 验证 。 网 格 无 关 性 验 
证 采用 的 参数 为 :H =2. 5c,y = 15°, J =0.3, Re = 
750。 为 了 避免 计算 域 尺 寸 对 计算 结果 的 影响 ,选择 
足够 大 的 流 场 尺 寸 , 流 场 的 外 边界 距 W2 BAAR 40c。 
在 表 1 中 给 出 采用 4 种 网 格 计算 后 的 W2 时 均 升 力 
系数 。 与 细 网 格 相 比 ,中 网 格 和 更 细 网 格 的 相对 误 
差分 别 为 1.95% 和 1. 12% 。 同 时 观察 图 4 可 知 , 升 
力 系数 和 推力 系数 随 网 格 单元 数 的 变化 并 不 明显 。 
因此 在 保证 计算 效率 和 满足 计算 精度 的 要 求 下 ,本 
研究 采用 细 网 格 进行 计算 。 图 5 中 给 出 了 计算 所 采 
用 的 网 格 示意 图 。 
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3 结果 与 讨论 


3.1 气动 力 分 析 


在 讨论 群 组 中 水 平 间距 互 对 扑 壳 群 组 气动 性 能 
影响 时 ,将 垂直 间距 y 固定 为 15" , 扑 翼 呈 斜 一 字 串 
联 排 布 。 图 6 绘制 了 群 组 中 各 个 扑 翼 的 时 均 气 动 性 
能 和 群 组 平均 气动 性 能 (mean ) 随 水 平 间距 的 变 
化 ,以 及 单 豆 ( 流 场 中 只 有 一 个 扑 辟 ) 的 气动 性 能 。 

首先 观察 群 组 中 个 体 的 气动 性 能 。W1 的 升力 
系数 在 所 有 水 平 间 距 (1.5c ~4c) PIa PAS, H. 
几乎 不 随 水 平 间距 H 的 变化 而 变化 。W2 和 W3 的 
升力 系数 呈现 先 增 后 减 小 的 趋势 。 当 及 =3c ~4c 范 
围 内 时 ,W2 和 W3 的 升力 系数 显著 高 于 单 汉 。 并 且 
观察 图 6(a) 和 图 6(c) ,可 以 发 现 群 组 中 个 体 的 升 
力 系 数 和 升力 效率 具有 相同 的 变化 趋势 。 

在 图 6(b) 和 图 6(d) 可 以 发 现 群 组 扑 恤 的 推力 
系数 与 推力 效率 具有 相似 的 变化 趋势 。W1 的 推力 
系数 和 效率 随 变化 较 小 ,并 与 单 辟 的 推力 系数 和 效 
率 相 近 。W2 和 W3 的 推力 系数 和 效率 随 着 水 平 间 
距 的 增加 而 明显 减 小 , 当 甩 >3. Oc 时 ,推力 性 能 小 于 


10 317 175 


(a) th 


0.659 2 


= 
“10 0.25 0.50 0.75 1.00 
©) z 
N (b) 
< 二 RZ ast E — aJ] 
Ò 4 RW 瞬时 气动 力 系数 
ig. nstantaneous aerodynamic coefficient 
Fig.4 I dy ffi 
ca of flapping wing W2 
© 表 1 网 格 无 关 性 验证 
CD Tab.1 Grid independence validation 
BERU 单元 数量 W2 (WC, 相对 误差 /% 
N: 3 673 357 0.6172 5.32 
2, 5 152 193 0.639 2 1.95 
<, 7 282 302 0.6519 


BAIMIINMMAS 


LO. SRL AY Po Hit ERTE E 


(b) 扑 性 表面 网 格 


Z 


Fig.5 Grid diagram 


5 网 格 示意 图 


PARE, 

从 扑 辟 群 组 的 整体 性 能 来 看 ,整体 的 平均 升力 
系数 和 效率 随 互 的 增加 呈现 先 增 大 后 减 小 的 趋势 ， 
平均 推力 系数 和 效率 随 着 有 的 增加 而 不 断 降低 。 并 
HERE H =3. 5c 时 整体 的 平均 升力 性 能 表现 最 
好 ,升力 系数 和 效率 分 别 比 单 豆 提高 62. 79% 和 
47.52% 。 此 时 在 推力 性 能 下 降 相 对 较 低 的 ,推力 系 
数 和 效率 分 别 下 降 了 13.35% 和 20. 43% if] H = 
1. 5c 时 ,整体 的 推力 性 能 最 好 。 与 单 缀 相 比 ,虽然 整 
体 的 平均 升力 效率 下 降 了 26. 56% ,但 整体 的 平均 
推力 系数 和 效率 分 别提 高 了 39.6% 和 5.54% ,升力 
系数 几乎 没有 下 降 。 这 些 结果 表明 , 群 组 中 的 个 体 
可 以 通过 改变 个 体 间 的 水 平 间距 来 获得 最 佳 升力 性 
能 或 推力 性 能 。 

由 以 上 分 析 可 知 ,水 平 间距 互 对 群 组 中 扑 翼 的 
气动 性 能 有 显著 影响 。 群 组 中 的 个 体 通 过 可 以 改变 
个 体 间距 ,达到 最 佳 升力 或 者 推力 性 能 。 这 可 能 是 
自然 界 中 马 类 在 迁徙 过 程 中 间 软 性 改变 其 运动 队 形 
MRAZ 。 


H 
(d) 


图 6 fe y =1S°RN EAA PPh Se EA EEEE 及 的 气动 性 能 


Fig.6 The aerodynamic performance of flapping wings at different horizontal distances H in a group with y = 15° 


3.2) 涡 系 结构 分 析 


Sey TBE AGE RAL PARAL 
扰 问 题 。 图 7 给 出 互 =1.5c( 最 佳 推力 性 能 ) #1 H = 
3 于 (最 佳 升力 性 能 ) 群 组 中 扑 翼 瞬 时 气动 力 曲线 。 
发 现 群 组 扑 必 升 力主 要 在 曲线 的 峰值 处 (7 =0.3 附 
近 7 和 低谷 (7 =0.7 附近 ) 处 产生 差异 ,而 推力 主要 在 
EAD r =0.7 时 刻 处 产生 差异 。 因 此 群 组 中 扑 翼 气动 
性 能 的 差异 主要 在 上 冲程 和 下 冲程 的 中 程 处 产生 。 
通过 研究 昆虫 飞行 发 现 ,稳定 的 前 缘 涡 (LEV ) 
可 以 增 大 昆虫 的 升力 *。 分 析 涡 的 脱落 机 理 和 变化 
规律 ,对 研究 群 组 中 存在 的 各 种 复杂 效应 及 其 个 体 
之 间 的 相互 作用 ,具有 重要 的 科学 价值 。 选 取 不 同 
水 平 间距 五 的 2 个 典型 时 刻 (升力 和 推力 最 大 时 
刻 ) 进 行 分 析 。 图 8 Ath TARH = 1. Se Al H = 
3.5¢ 三 种 情况 下 ,r =0.27 和 7=0.67 时 0.5R 处 截 
面 涡 量 图 。 首 先 ,观察 下 冲程 (图 8,7 =0.27)。 在 
7=0.27 时 ,比较 了 五 =1.5c o BARS Ath ee 
群 运动 时 的 涡 量 图 ,发 现 Wl 的 后 缘 涡 (TEV) 非常 
接近 W2 的 前 缘 , 使 W2 AY LEV EAER, FAEH 
于 W2 的 TEV 的 诱导 , W3 的 上 表面 也 形成 了 更 大 
HJ LEV, XF H =3. 5c 的 情况 (图 8(c)), 在 W2 和 


W3 的 上 表面 也 观察 到 了 更 大 的 LEV, (A H Pe 
的 分 离 TEV 到 达 后 翼 并 与 后 经 相互 作用 需要 更 多 
的 时 间 , 因 此 万 =3. 5c 情况 下 W2 和 W3 的 正 升力 
峰值 出 现 的 时 间 比 万 =1.5c 情况 要 晚 ( 图 7)。 在 
Boschisch' "等 人 的 实验 研究 中 也 观察 到 了 随 着 水 
平 间 距 的 增加 相 邻 扑 玻 之 间 的 延迟 相互 作用 。 

在 上 冲程 时 ,在 下 冲程 中 扑 缀 产生 的 TEV 已 从 
后 缘分 离 并 运动 到 下 游 ( 图 8(a),7 =0.67)。 在 
H=1.5c 中 ,W1 分 离 的 TEV 被 W2 所 捕获 ,从 而 加 
强 了 W2 上 CCW LEV 的 形成 ,这 导致 W2 的 正 推力 
和 人 负 升 力 的 增 大 ( 见 图 7(b))。 同 样 , W2 分 离 的 
TEV 也 被 W3 捕获 ( 见 图 8(p),r=0.67) ,增强 了 
W3 的 CCW LEV( 注 意 :W1 分 离 的 TEV 和 W2 的 
LEV 沿 同一 方向 旋转 ( 逆 时 针 , CCW))。 在 二 = 
3. Se 中 ,由 于 扑 避 的 水 平 间 距 相 对 较 大 ,前 豆 分 离 的 
TEV 迅速 消散 ,W2 和 W3 SAJATI TEV 的 影响 较 
小 ,解释 了 五 =3.5c 时 群 组 中 扑 翼 的 气动 力 与 上 剖 
程 时 单 辟 的 气动 力 为 何 相似 (图 7(c) 和 图 7(d) )。 
观察 图 8 中 7 =0.67 时 刻 ,发 现 群 组 中 的 后 嗓 由 于 
ALF A BEA SAY TEV 影响 ,W2 和 W3 的 CCW LEV 
Ve PE ELT. E W2 和 W3 AY Spe bE eT EY 
TEV 完整 性 遭 到 破坏 ,因此 在 r=0.67 中 W2 All W3 
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数 可 以 看 出 ,推力 产生 差异 的 时 刻 主 
6 


0.75 


0.50 
T 


(c) H=3.5c 


0.75 


ChinaXivA ERAT! 
中 程 处 ,所 以 五 =3.5c 时 W2 All W3 的 推力 系数 小 


0.50 
T 


(b) H=1.5c 


0.75 1.00 


0.25 


0.50 
T 


(d) H=3.5c 


0.75 


图 7 


在 不 同 水 平 间 距 有 下 , 群 组 中 每 个 扑 翼 和 和 


权 瞬 时 升力 和 推力 系数 曲线 


Fig.7 Comparison of the lift and thrust coefficient over a flapping cycle for each individual wing and 
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the single wing at different horizontal spacing H 


ee roe so 较 大 的 水 平 间距 下 , RSE AY DLA TL, A 
500 捉 没 有 发 生 。 
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